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Abstract 

We define the entropy S and uncertainty function of a squeezed system interacting 
with a thermal bath, and study how they change in time by following the evolution 
of the reduced density matrix in the influence functional formalism. As examples, we 
calculate the entropy of two exactly solvable squeezed systems: an inverted harmonic 
oscillator and a scalar field mode evolving in an inflationary universe. For the inverted 
oscillator with weak coupling to the bath, at both high and low temperatures, S — > 
r, where r is the squeeze parameter. In the de Sitter case, at high temperatures, 
S ^ {1 — c)r where c = jq/H, 70 being the coupling to the bath and H the Hubble 
constant. These three cases confirm previous results based on more ad hoc prescriptions 
for calculating entropy. But at low temperatures, the de Sitter entropy S — > (1/2 — c)r 
is noticeably different. This result, obtained from a more rigorous approach, shows 
that factors usually ignored by the conventional approaches, i.e., the nature of the 
environment and the coupling strength betwen the system and the environment, are 
important. 
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1 Introduction 



In discussing the conceptual problems of entropy generation from cosmological particle cre- 
ation Ijl], one of us was confronted in the early 80's by the following apparent paradox: 
on the one hand common sense suggests that entropy (S) is given by the number (A^) of 
particles produced {S ~ A^^ for photons). On the other hand, theoretically, for a free 
field, particle pairs created in the vacuum will remain in a pure state and there should 
be no entropy generation. Inquiry into this paradox led to serious subsequent investiga- 
tions into the statistical properties of particles and fields. In 1984, Hu pointed out that 
the usual simplistic identification of entropy with the number of particles present is valid 
only in the thermodynamic-hydrodynamic regime, where interaction among particles and 
coarse-graining can lead to entropy generation. This aspect was discussed later by Hu and 
Kandrup using a statistical mechanics subdynamics analysis. The more intriguing case of 
entropy generation for free fields was addressed by Hu and Pavon 0. They suggested that 
an intrinsic entropy of a (free) quantum field can be measured by the particle number (in 
a Fock space representation) or by the variance (in the coherent state representation). The 
entropy of a (free) quantum field is non-zero only if some information of the field is lost or 
excluded from consideration, either by choosing some special initial state and/or introduc- 
ing some measure of coarse- graining. For example, the predicted monotonic increase in the 
spontaneous creation of bosons is a consequence of adopting the Fock space representation 
which amounts to a random phase initial condition implicitly assumed in most discussions of 
vacuum particle creation. (The difference of spontaneous and stimulated creation of bosons 
versus fermions was first pointed out by Parker and discussed in squeezed state language 
by Hu, Kang and Matacz 0). The relation of random phase and particle creation was 
further elaborated by Kandrup [§. 

Following these early discussions of the theoretical meaning of entropy of quantum fields, a 
recent surge of interest on this issue was stimulated by the work of Brandenberger, Mukhanov 
and Prokopec (BMP) 0, Gasperini and Giovannini (GG) |10| and others on the entropy 
content of primordial gravitons. The language of squeezed states for the description of 
cosmological particle creation was introduced by Grishchuk and Sidorov [|ll|. Though the 



physics is the same ||^, |T^] as originally described by Parker [|T| and Zeldovich |Q, the lan- 
guage brings closer the comparison with similar problems in quantum optics, which shares 
many interesting theoretical and practical issues |jl3|. BMP suggested a coarse-graining of 
the field by integrating out the rotation angles in the probability functional, while GG con- 
sidered a squeezed vacuum in terms of new variables which give the maximum and minimum 
fluctuations, and suggested a coarse-graining by neglecting information about the subfluc- 
tuant variable. Keski-Vakkuri studied entropy generation from particle creation with many 
particle mixed initial states Matacz [|1^] considered a squeezed vacuum of a harmonic 
oscillator system with time-dependent frequency, and, motivated by the special role of co- 
herent states, modeled the effect of the environment by decohering the squeezed vacuum 
in the coherent state representation. Kruczenski, Oxman and Zaldarriaga |16| also used a 
procedure of setting off-diagonal elements in the density matrix to zero before calculating 
the entropy. Despite the variety of coarse-graining measures used, in the large squeezing 
limit (late times) these approaches all give an entropy of 5 = 2r per mode, where r is the 
squeezing parameter. This result which gives the number of particles created at late times 
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agrees with that obtained in the original work of Hu and Pavon . 

Noteworthy in this group of work is that the representation of the state of the quantum 
field and the coarse-graining in the field are stipulated, not derived. What is implicitly 
assumed or grossed over in these approaches is the important process of decoherence - the 
diminution of the off-diagonal components of a reduced density matrix in a certain basis. It is 



a necessary condition for realizing the quantum to classical transition |jT^. The deeper issues 
are to show explicitly how entropy of particle creation depends on the choice of specific initial 
state and/or particular ways of coarse-graining, and to understand how natural or plausible 
these choices of the initial state representation or the coarse-graining measure are in different 
realistic physical conditions [1^ .tl To answer these questions, one needs to work with a more 
basic theoretical framework, that of statistical mechanics of quantum fields. In recent years 
we have approached the decoherence and entropy / uncertainty issues with the quantum open 



system concept and the influence functional formalism [^. The purpose of this paper 



is to study the entropy and uncertainty of quantum fields using the statistical mechanics of 
squeezed quantum open systems as illustrated by quantum Brownian motion models. 

In the quantum Brownian motion paradigmic depiction of quantum field theory studied 
in the series of papers by Hu, Paz, Zhang [^, ^] and Hu and Matacz [Q, the system 



represented by the Brownian particle can act as a detector (as in the influence functional 



derivation of Unruh and Hawking radiation |25, p6[), a particular mode of a quantum field 



(such as the homogeneous infiaton field), or the scale factor of the background spacetime (as 
in minisuperspace quantum cosmology), while the bath could be a set of coupled oscillators, 
a quantum field, or just the high frequency sector of the field, as in stochastic inflation. The 
statistical properties of the system are depicted by the reduced density matrix (rdm) formed 
by integrating out the details of the bath. One can use the rdm or the associated Wigner 
function to calculate the statistical average of physical observables of the system, such as 
the uncertainty or the entropy functions. The von Neumann entropy of an open system is 
then 

5* = -tr Pred In Pr-ed (1-1) 

The uncertainty function measures the effects of vacuum and thermal fluctuations in the 



environment (at zero and finite temperature) on the observables of the system [27, 28]. The 
increase of their variances due to these fluctuations gives rise to the uncertainty and entropy 
increase. The time-dependence of the uncertainty function of an open system measures the 
varying relative importance of thermal and vacuum fluctuations and their roles in bringing 



about the decoherence of the system and the emergence of classical behavior [p7| , p8 |. 

The entropy function constructed from the reduced density matrix (or the Wigner func- 
tion) of a particular state measures the information loss of the system in that state to the 
environment (or, in the phraseology of [^], the 'stability' characterized by the loss of pre- 
dictive power relative to the classical description). One can study the entropy increase for 
a specific state, or compare the entropy at each time for a variety of states characterized by 
the squeeze parameter. The time scale of entropy increase, when entropy arises from particle 



^This includes conditions when, for example, the quantum field is at a finite temperature or is in dise- 
quilibrium, interacting with other fields, or that its vacuum state is dictated by some natural choice, e.g., 
in the earlier quantum cosmology regime such as the Hartle-Hawking boundary condition leading to the 
Bunch-Davies vacuum in de Sitter spacetime. 
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creation from the vacuum, should be comparable to the decoherence time, which, for a high 
temperature bath, is very short. Interaction with the environment also changes its dynamics 
from strictly unitary to dissipative, the energy loss being measured by the viscosity function, 
which governs the relaxation of the system into equilibrium with the environment. The 
entropy function for such open systems can also be used |2^, |29[ as a measure of how close 
different quantum states can lead to a classical dynamics. For example, the coherent state 
being the state of minimal uncertainty has the smallest entropy function and a squeezed 
state in general has a greater uncertainty function [27]. One can thus use the uncertainty to 



measure how classical or 'nonclassical' a quantum state is. 

Using this first-principle approach for the calculation of the entropy function leads to 
more reliable results. With regard to the issue of entropy of quantum fields raised at the 
beginning, we can now ask, what is the difference of our more vigorous definition and that 
defined earlier with more ad hoc prescriptions? 

Foremost, the differences in design are obvious: the entropy of M, 0, 110] and others refers 



to that of the field, and is obtained by coarse-graining some information of the field itself, 
such as making a random phase approximation, adopting the number basis, or integrating 
over the rotation angles. The entropy of [^, ^ refers to that of the open system and 



is obtained by coarse-graining the environment. Why is it that for certain generic models in 
some common limit (late time, high squeezing), both groups of work obtain the same result? 
Under what conditions would they differ? Understanding this relation could provide a more 
solid theoretical foundation for the intuitively-argued definitions of field entropy. 

At the formal level, supposing we have some system which has been decomposed into two 
subsystems, it can be shown [PU| that between the entropies 81,82 of the two subsystems. 



and that of the total system, 5*12, a triangle inequality holds: 

\8i-82\<8u< 81 + 82 (1.2) 

In particular, if the total system is closed and so in a pure state, then it has zero entropy, 
so that the two subsystems necessarily have equal entropies.! Hence, asking for the entropy 
change of a system is equivalent to asking for the entropy change of the environment it 
couples to, if the overall closed system is in a pure state. Now consider the case of the 
system as a detector (or a single mode of a field) and the environment as the field. The 
information lost in coarse-graining the field which was used to define the field entropy in 
the above examples is precisely the information lost as registered in the particle detector, 
which shows up in the calculaton of entropy from the reduced density matrix. The bilinear 
coupling between the system and the bath as used in the simple quantum Brownian motion 
models also ensures that the information registered in both sectors are directly commutable. 
This explains the commonalities. However, not all coarse-graining and coupling will lead to 
the same results, as we shall explicitly demonstrate in some examples. 

Another important feature of the entropy function obtained in our present investigation 
which is not at all clear in earlier studies is that it depends nonlocally on the entire history 



^This could be the reason why the derivation of black hole entropy (see the recent review of Bekenstein |3T) ) 
can be obtained equivalently by computing the entropy of the radiation (e.g., p2[) emitted by the black hole, 
or by counting the internal states (if one knows how!) of the black hole (e.g., [|33[). Physically one can view 
what happens to the particle as a probe into the state of the field. The application of open-system concepts 
to black hole entropy is a very fruitful avenue p4| . 
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of the squeezing parameter. This can be seen from the fact that the rate of particle creation 
varies in time and its effect is history dependent Existing methods of calculating 

the entropy generation give results which only depend on the squeezing parameter at the 
time when a particular coarse-graining (or dropping the off-diagonal components of the 
density matrix) is implemented. These ad hoc choices (of coarse-graining and the time it is 
introduced) affect the generality of the earlier results. 

The plan of this paper is as follows. In Sec. 2 we give a brief summary of a squeezed 
quantum system, using a general oscillator Hamiltonian as an example. The notation is that 
of 0, |2^ . In Sec. 3 we give a brief summary of open quantum systems in terms of influence 



functionals , following the treatment of |2^ . Readers familiar with these background 



material can go directly to Sec. 4, which contains the central material for the derivation of 
entropy and uncertainty functions. In Sec. 5 we apply these formulas to an oscillator system, 



recovering en route the earlier results of |27, Bq] for uncertainty at finite temperature, and 



of on entropy of coherent states. In Sec. 6 we apply our result to the consideration of a 
scalar field in a de Sitter universe. We show the conditions where one recovers the S = 2r 
result of all previous work, and more significantly, the cases when they differ. We give a 
short discussion of our findings in Sec. 7. The Appendices contain details of derivations. 

2 Squeezed Systems 

2.1 Squeezed states and density matrices 

Consider the general oscillator Hamiltonian 

2 

H{t) = fit) J + f*{t)^ + h{t){a^a + 1/2) + d{t)a + d*{t)a) + g{t) (2.1) 

where c?, /, h are arbitrary functions of time. The propagator for this has been calculated 
in and is 

f/(t, t,) = 5(r, c^)R{e)D{p)e'"-\^\'" (2.2) 
where p, w are defined in terms of the coefficients appearing in if, and 

D{p) = exp{—p*a — h.c.) 
R{9) = exp-ie{a^a + l/2) 
S{r,(f)) = exp(re-2^V/2 - h.c.) (2.3) 

are the displacement, rotation and squeeze operators respectively. Suppose we start with 
a simple harmonic oscillator with lagrangian 

L = ^ _ n^x^) (2.4) 

If we construct a gaussian state in the position basis, with initially the same width ctq as 
that of the ground state of such an oscillator, displaced by some arbitrary amount and with 
a phase proportional to x, we find this to be an eigenstate of the lowering operator, and is 
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called a coherent state. Suppose we locate the point ((x), (p)) in phase space and draw an 
ellipse about this point, the lengths of whose axes being the uncertainties Ax^, Ap^. Then as 
the oscillator evolves this uncertainty ellipse revolves about the origin with angular speed Q. 

A squeezed state is again such a state, but with an arbitrary initial width a. We find 
that as the oscillator evolves the uncertainty ellipse again revolves about the origin, but its 
axes change length and it can also rotate about its own centre. 

It turns out that the squeeze parameter r is related to the width of such a state: 



In 



2Mn 



(2.5) 



Hence a coherent state has r = 0, or zero squeezing. A gaussian that initially has a width 
smaller than ctq will evolve to a squeezed state with some r > 0. We can generate a squeezed 
state by applying S{r, (p) to the ground state of the simple oscillator. Consider the new 
operator 

b = U^aU = aa + P*a^ (2.6) 

where it turns out that 



a 



-ie 



e chr 



(2.7) 



Going from a to 6 is then just a Bogoliubov transformation, and so a, (3 become Bogoliubov 
coefficients for our system. Their equations of motion are 

a = —iha — if* [3 

(3 = ifa + ihp (2.8) 
a{ti) = 1 , Piti) = 

where /, h as defined in the hamiltonian (|2.1|) are calculated from the general system la- 
grangian. This lagrangian has time dependent mass and frequency, and we will also allow it 
to have a time dependent cross term denoted 2S{t): 



(x^ + 2£{t)xx - n^{t)x^) 



(2.9) 



Then /, h are given by |24 



/ 
h 



K 



K 


+ 2iS 


M 






M. 





(2.10) 



and K is an arbitrary positive constant that can be chosen to simplify the relevant equations. 

In the next section we shall find that the quantity of much importance to our work turns 
out to be the sum of the Bogoliubov coefficients, X = a + p. It follows from ( ^.8] ) that X 
satisfies the classical equation of motion for the system: 



Q' + S + 



X = 



(2.11) 
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with initial conditions 

X{U) = 1 ; X{U) = ^-S{U) (2.12) 

With this result, the usual task of finding the Bogoliubov coefficients a, (3 from two coupled 
first order differential equations is reduced to that of solving one second order equation for X. 

2.2 Squeezing an inverted harmonic oscillator 

For an inverted oscillator, i.e. one with < 0, at late times r is expected to blow up. In 
that case we can calculate it from (U-7\) as follows, h 



\a\ \P\ e72 (2.13) 

so that 

r^ln(2|a|) (2.14) 
Rather than use (|2.8|) to calculate a, once we have X we can extract a from it. This is done 



by writing, from (|2!8 



X = a + p 

X = i{f -h)a + t{h- r)f3 (2.15) 
and solving for using (|2.10|) : 

We can follow the behaviour of r, 0, 6 by writing (|2.8|) in terms of the squeeze parameter, 
with / = \f\e'': 

r = \f\sm{2(f) + e) 

= |/|coth2rcos(20 + e) 

e = /i- l/l thrcos(20 + e) (2.17) 

These equations are useful for numerical work. They also tell us of the existence of constant, 
and so possibly attractor, solutions for 0, 6. If we set r — oo then the equations for 0, 6 
become 

^ = -0 = /i - |/|cos(20 + e) (2.18) 

1. Suppose there exist some 9 and such that 9 = (p = 0. Then h = \ f\ cos(20 + e), so 
that \h\ < \ f\. Thus, since h is real, we have h"^ < and from ( |2.10| ) this inequality 
is true if and only if < 0. 

2. Conversely suppose < 0. Then by the previous argument, \h\ < |/|, or —1 < 
h/\f\ < 1. Thus there must exist some such that cos(20 + e) = h/\f\. From ( p.l8| ) 
we see that for this value of 0, 9 = 6 = 0. 
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In other words, there will exist constant solutions for 0, 9 if and only if < (the oscillator 
is "inverted"). Of course, this doesn't reveal whether these constant solutions are attractors. 
Numerically solving (|2.17 ) with fi^ < 0, for various VL and k, shows that 0, 9 apparently 
do always quickly tend toward constants, always accompanied by one of r — >■ ±oo. 

We note that it's common to eliminate the cross term in the action by adding a surface 
term: 



L 



— (x^ + 2Sxx - - --{MSx') 

2 V ^ 2dV ' 



M 

T 



x^-{a' + 



MS 



S\x' 



(2.19) 



Although this leaves the classical equation of motion unchanged, it will change the squeeze 
parameters. In this paper we leave the cross term in our lagrangians. 



3 Open Systems 



3.1 Influence functional theory 



20 



as 



The influence functional (IF) formalism was first introduced by Feynman and Vernon 
a way of deducing the infiuence of an environment on some system of interest. It was later 
applied by Caldeira and Leggett ||2^ to the high temperature limit of a model where both 
system and environment are composed of static oscillators, that is, having time independent 
frequency. A comprehensive review is given by Grabert et al (in pO[|). 

In these earlier works, the infiuence functional for quantum Brownian motion has only 
been derived for Markovian processes corresponding to coupling to a high temperature ohmic 
bath. An exact master equation for non-Markovian processes is recently derived by Hu Paz 

mm 



and Zhang |2^, |23| (see also 0)- Hu and Matacz ||2^ obtained the master equation for 
system and bath oscillators with time- dependent frequencies, a result readily generalizable 
to quantum fields. Stochastic properties of interacting quantum field theory are discussed in 
pO| , 41 1 . Most work in this area since Feynman and Vernon has assumed a bilinear system- 
bath coupling, which yields an exact analytic form for the infiuence functional. Recently, 
weak nonlinear couplings p3[ have also been considered using perturbation theory borrowed 



from field theory. 

The language of infiuence functionals was developed in the context of non-equilibrium 
statistical mechanics, but can be generalized to field theory (see e.g., |41]). In fact it can 



be shown to be formally equivalent to the Schwinger-Keldysh closed time path (CTP) 
Stochastic field theory based on the IF and CTP has since been applied to 



formalism [21 



semiclassical gravity |Q and infiationary cosmology problems 



In this paper we further develop the work of p4| by considering a squeezed system 
coupled bilinearly to a static bath (oscillators with time-independent frequencies), but with 
a time-dependent coupling constant. We also lay out the groundwork for calculating such 
quantities as entropy and uncertainty as well as fiuctuations and coherence, for the purpose 
of this paper, and a later one on the de Sitter universe p5|. 
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3.2 Propagator for the density matrix 

The primary object we wish to consider is the evolution of the reduced density matrix of our 
system via the Feynman- Vernon influence functional method. This has been discussed at 



length in [Q; we describe it here in order to establish the notation, and just state its main 
results without deriving them. 

Again consider our system described by x which interacts with its environment q through 
some interaction. The combined action is 

S[x, q] = S[x] + SE[q] + Sint[x, q] (3.1) 

We require the reduced density matrix of the system at time t. This is found by tracing out 
the environment: _ 

Pr{xx't)= / dqp{xqx'qt) (3.2) 



The full density matrix p{xqx' qt) evolves unitarily. Suppose we expand it using complete- 
ness relations and then path integrals: 

p{xqx' qt) = {x qt\p\x' qt) 

dxidqi J dx[dq[ {x qt\xiqiQ){xiqiQ\p\x[q[Q){x[q[Q\x' qt) 

dxidqi f dx[dq'i r Dx T e^^[^''?lp(a;i g- 0) T Dx' f Dq' e-^^l^''-^'^ 

J Jxi Jqi Jx'^ Jq'. 

= J dxi dqi J dx[ dq[ J{x qx' qt\xiqix[ g- 0) p{xi qi x[ q[ 0) (3.3) 

where J is seen to be an evolution operator for the entire system plus bath. Now to allow fur- 
ther calculation we make the assumption that the system and bath are initially uncorrelated, 
i.e. 

p{xi qi x'i q[ 0) = psys{xi x- 0) pE{qi q'i 0) (3.4) 

(Initial conditions with correlations have also been considered by |^). In this case we are 
able to rearrange the order of integration to write the reduced density matrix in the following 
way: 

j-^iY" X t ^ — / dx 2 dx ^ J p ^ X X 1 1 X ^ X ^ ^ 1^ gy g ^ X i X ^0^ ^3.5^ 



where the evolution operator for the reduced density matrix is defined by 

Jr{xx't\xix[0) = rDx rDx'e^^[^l-^^[^'l F[x,x'] (3.6) 

Jxi J x'- 

and F[x, x'] is the so-called influence functional: 

F[x, x'] = Jdq dqi dq[ psiqi q'i 0) J' Dq e*^^=M+iS»t[x,q] j"^ j^^' ^~iSE[g']-iS„uW,g'] (3 7) 

We can also write the influence functional in a basis-independent form as follows. First we 
write the path integrals as propagators 

F[x,x'] = [dq dq, dq[ pE{qii^) {q\U{t)\q,) {q[\U'\t)\q) (3.8) 



where U{t), U'{t) are the propagators for Se[(1\ + Sint[x, q\ and Se[(1\ + Sint[x' , q\ respectively. 
Then upon integrating over g, qi and writing the remaining integral as a trace, we obtain: 



F[x,x'] = tr U{t) pE{0)U'\t) 



(3.9) 



Using this form to calculate the influence functional was done earlier in ||2^. Here we just 
list the result: if we use sum and difference coordinates defined by 



J: = {x + x')/2 , A = x-x' 



(3.10) 



then the influence functional can be written in terms of two new quantities, the "dissipation'' 
fi{s,s') and "noise" i/(s,s'): 



(3.11) 



F[x,x'] = exp —— ds ds' A{s) z/(s, s') A(s') + 2/i(s, s') 2E(s') 
n Jo Jo L 



Thus the influence of the environment is completely invested in the dissipation and noise. 



3.3 Evolution of the reduced density matrix 

Suppose now that we work within the context of quantum brownian motion, using the 



notation of That is, our system is modeled by an oscillator with time dependent mass, 
cross term and natural frequency. This interacts bilinearly with an environment modeled in 
the same way, the total lagrangian being 

S[x, q] = S[x] + SeIq] + Sint[x, q] 



J' ds I (x' + 2S{s)xx - n^{s)x' 



+ E 



m„(s) 



(ql + 2en{s)qnqn - ^n(s)?r 



+ ^[-c(s)xg„] (3.12) 



where the particle and the bath oscillators have coordinates x and g„ respectively. 

We wish to start with some initial system density matrix psys {xi x\ 0) and evolve it us- 



ing (|3.5| ). As described in p4|, is calculated using the standard path integral approach. 
Using the sum and difference coordinates deflned in ( |3.10D , the classical paths followed by 
the system, E^;, A^;, can be written in terms of more elementary functions u^v: 



Then it can be shown that the superpropagator J,, is equal to 

I&2 



(3.13) 



27ih 



exp 



h 



(6iSA - 62SA, + fcsS.A - 64S.A, 



- -(anA^ + aisAiA + aasA^ 



(3.14) 
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The functions bi 64 can be expressed as 



bi{t,ti) = M{t)u2{t) + M{t)£{t) 

b2it,ti) = M{U)u2{ti) 

b3{t,ti) = M{t)u,{t) 

b^{t,U) = M{U)u^{ti) + M{ti)£{U) (3.15) 



while the functions defined by 

aijit.U) = , , . ds ds' Vi{s) z/(s, s') Vj{s') (3.16) 
1 + Oij Ju Jti 

The functions Ui — > V2 are solutions to the following equations (dropping subscripts on u, v): 
M ( ■ M \ 2 



Hs) + -v+io^+8 + —S\v- ^ ds' /i(s, s') v{s') = (3.18) 
subject to the boundary conditions 

uiiti) = vi{ti) = 1 , ui{t) = vi(t) = 

U2iU) = V2{U) = , M2(t) = V2it) = 1 (3.19) 

3.4 Propagator for the reduced density matrix: ohmic envi- 
ronment 

To proceed further we need explicit expressions for an — > 64. These are expressed in terms 
of Ml — > f2, which in turn come from solving ( p.l7| , p.l^ ). To solve these equations we need 
to know the dissipation fi of the environment. 



M{s) Ju 
2 

M{s) 



The noise and dissipation can be calculated from p^, eqns 2.18, 2.19]. We choose the 



bath oscillators to be simple harmonic, that is, static with no cross term, since this turns 
out to correspond to the simplest form of dissipation: local. For such an environment the 
dissipation and noise can be shown to be 

POO 

fi{s,s') = / du I{u,s,s')lm[X{s)X*{s')] 
Jo 

v{s,s') = / du I{uj,s,s')coth—Re[X{s)X*{s')] (3.20) 
JO 2T 

where by T we will always mean ksT /h\ X is the sum of the Bogoliubov coefficients for the 
bath oscillators and / is the "spectral density" , a function defined by 

/(^, s') = ^^4^ E ^ - ^n) (3.21) 
2k 
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which encodes information of the action of the environment on the system. In general the 
spectral density can be described by some function of , where j is set by the particular 
environment being modeled. The case of j = 1, a so-called "ohmic" environment, is a 
borderline between the super-ohmic case (j > 1) — which models weak damping — and the 
subohmic case (j < 1) modelling strong damping. We can in effect consider both damping 
extremes by taking an ohmic environment together with some strength 70 which can be 
altered from zero, for a free system, up to higher values. 

Also, by considering the continuum limit of the coupling constant, it can be shown that 
this constant's independence of n also leads to an ohmic environment; so we will only consider 
spectral densities of the following form: 

I{uj,s,s') = —ujc{s)c{s') (3.22) 

TT 

For a general lagrangian the sum of the Bogoliubov coefficients X will be complicated; 
however we have simplified our calculations by taking the bath to be composed of unsqueezed 
(i.e. coherent) static oscillators with unit mass. For this type of bath the dissipation and 
noise can be calculated for an arbitrary bath temperature; we use the integral form of the 
noise as being easier to work with: 

fi{s,s') = 2-foc{s)c{s')6'{s-s') 
270 



i/(s, s' 



c{s)c{s') [ (jj coth-^ COS (jj( s — s') d(jj (3.23) 
71 Jo 2T 



In the high temperature limit the noise becomes white, that is it tends toward a delta 
function. 



4 Entropy and uncertainty, fluctuations and coher- 
ence 

4.1 Initial and final states 

Assume the systems are initially in the vacuum state, so that their density matrix is gaussian. 
So we start with an arbitrary gaussian reduced density matrix 

Pr{xi x[ U) oc e-«"?+^"-"--«*<' (4.1) 

and propagate it by using ( p.5| , p.l4| ) to give 

-AA2-2jBAS-4CS2 



p,(x, x', t) = iVe-^^^-2*^^^-4^'^^ (4.2) 



where we have used the same A, B and C notation of ||T5[, and with ^r^ii the real and 
imaginary parts of ^: 



N = 2^C/^ 

A = 022 + {[(2^,. + x)/4 + an] bl + (2^^ + 64) h - (2^. - x)al2} 
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1 



D = 4|er-x' + 4(2e.-x)aii+46&4 + &4 (4-3) 

These expressions form the basis of our later calculations. The quantity we are focusing on 
is the reduced density matrix, ([4.2|) , using the expressions in ( [4.3|) . These in turn use ( |A.14|) , 
which depends on our obtaining X, the sum of the Bogoliubov coefficients for the effective 
oscillator. 



4.2 Entropy from the reduced density matrix 

The entropy of a field mode has been calculated by Joos and Zeh [|^. It can be derived 
from the reduced density matrix at time t by using ( |1 . 1| ) , and is given by 



S = — [wlnw + (1 — w) ln(l — w)] ~ 1 — Inw iiw 



where 



w 



2JC/A 



1 + ^C/A 

The linear entropy is often more useful to work with owing to its simplicity: 

Sun = -trp2 = -JcJa 



and S = ^ oo is equivalent to Sun = 
we have 

5 ^ - In \Siin\ + 1 - In 2 



1^0, both strictly increasing. Then if S, 



lin 



i.e. S, 



lin 



—e 



i-s 



/2 



(4.4) 
(4.5) 

(4.6) 


(4.7) 



As an example, suppose we have a system in an initially pure gaussian state (x = 0), so that 
noise and dissipation are absent: 70 = 0. In this case, from ( |3.23| , |A.14| ) we have 



so that 



an = CL12 = 022 = 
gives C/A = 1 and hence from (^^) = as expected. 



(4i 



4.3 Fluctuations and coherence 

A clearer picture of the dynamics of a closed and open system can be obtained if we rotate 
the phase space axes so that the density matrix can be expressed in terms of the so called 
super- and subfiuctuant variables. (Alternatively, we are rotating the Wigner function in 
phase space so as to eliminate the cross term there). Call these variables u,v, expressed as 
real linear combinations of q,p [they have nothing to do with the u,v of ( p. IS] )]. We fix the 
linear combinations such that one variable {u, the superfiuctuant) grows exponentially while 
the other decays exponentially. In the case of no coupling to the environment we proceed 
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by expressing (m^), (f ^) in terms of (g^), {qp + pg), (p^), and then substituting for these the 
standard squeezed state resuhs [l^. This enables us to write 



ne 



2r 



-2r 



2k 



(4.9) 



These relations fix u, v in terms of q,p, and we now use the same transformation for the case 
of nonzero dissipation: 



u = — Ksin0g + COS0P 

, sin0 
V = cos (pq -\ p 

K 



What we wish to do is take a density matrix in position, 
Consider first of all calculating p{u,u'): 



(4.10) 

and write it in the m, v basis. 



p{u,u') = / {u\q) p{q,q') {q'\u) dqdq' 



(4.11) 



We need {u\q). This can be found by solving the p.d.e which follows by quantising ( [4.10 ) 
and applying both sides to {q\u)\ 



u{q\u) = {—Ksincf) q — icoscf) dg) {q\u) 



which has solution 



{q\u) = /(M)exp- 



COS( 



K sm (p q 



+ qu 



(4.12) 



(4.13) 



for some function f{u) to be determined [unrelated to ( p.lO| )]. We determine f{u) by redoing 
this calculation with the roles of q and u interchanged; since [f , u] = i, we have 



q{u\q) 



sm u 



K 



+ icos(j)du {u\q) 



Solving this determines f{u) and allows us to finally write (up to a phase 

1 i 



{q\u) 



exp ■ 



A/27rcos0 cos 



Similarly we find 



{q\v) 



2tx sin I 



exp 



IK 



sm I 



K sin q^ sin v? 



— cos (f) q^ cos (f) v"^ 
V qv 



(4.14) 



(4.15) 



(4.16) 



Now, suppose we start with a gaussian density matrix as in ( |4.2| ). We can then easily change 
bases using ( ^TTj ^l5| , ^l6D to get, with 



7 = - cot > 



sm 



AAC +{B- -if 



(4.17) 
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A = 



4AC + (47a + B--ff 
4^^ 



(4.1^ 



p{v,v') 



C 
TxaX 

C 



exp 




4aA 

-1 



AAl + 2«(47a + B- 7)A„E„ + AC^l 



exp — 
TTcr 4cr 



AA^ - 2z(47(T + 5 - 7)A„S„ + 4CS 



(4.19) 



where we have used sum and difference variables, e.g. £„ = (« + u')/2, A^ = u ~ u', and 7 
has no relation to 70. 

We can show that in the absence of a bath, these matrices reduce to the expected ones for 
a squeezed vacuum. First, in the g-representation the density matrix of a squeezed vacuum 
is known to be |57] 



piq,q') oc 



-K 



1 + e^^'^tanhi 



(4.20) 



2 1 - e2*<^ tanh r 

If we write p{q, q') in terms of sum and difference coordinates and compare with the defini- 



tions of A, B, C in (^), we find 

A = C 

B 



1 — tanh r 



Substituting these into ( |4.19| ) gives 

p{u,u') -- 
p{v,v') -- 



4 1 — 2 cos 20 tanh r + tanh^ r 

K sin 20 tanh r 
1 — 2 cos 20 tanh r + tanh^ r 



e — e 

: exp 



(4.21) 



K —KC 

—e exp 

TT 



-2r 

— (W^ + «'^) 

2r 

— (v"^ + V 



(4.22) 



These are the expected results, as can be seen by the fact that with p, q replaced by u, v 
respectively, they are produced when is set to zero in p{p,p') and p{q,q'). 

Measures of fluctuations and coherence 

Returning to the general case of dissipation, the fluctuations in u and v are calculated from 
the density matrices: 



{u^) — (u)^ = / u^p{u,u) du — 



u p{u, u) du 



2C 



A 2 cr 



(4.23) 



and both of these are just equal to 1/2 divided by the coefficient of —T? in their density 
matrix. 
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As a measure of coherence we note that a large coefficient of — means that the density 
matrix is strongly peaked along its diagonal, i.e. there is very little coherence in the system. 
A measure of coherence was defined in |48| as a squared coherence length L^, equal to 1/8 
divided by the coefficient of — A^, so that a large means a high degree of coherence in the 
system. With this definition of L^, ( |4.19| ) gives 



^ r2 

2 A' " 



a 

2A 



(4.24) 



We can also relate the coherence lengths and fluctuations to the entropy of the system (see 
section [4.2| for definitions) . We can write 

(A note of caution: linear entropy is negative by definition in order for it to increase with 
S. Then as Sun increases, will decrease). Also the uncertainty relation for m, v becomes, 
from (|4T7|, CT, ra) : 



An^At;^ 



C2 



1 (47cr + 5-7)^ 
4 ^ IQAC 



For the free field the last term in the square brackets is zero while Sun 
so that Am Aw = 1/2. 



(4.26) 

— 1 (since S* = 0), 



5 Entropy and uncertainty of oscillator system 

We can now demonstrate how the previous results are used. In the simplest cases, such as a 
static oscillator coupled to a thermal bath of static oscillators, with a static ohmic coupling, 
the entropy is easily compared with known results in equilibrium statistical mechanics. From 
section |3^ , we know that this case has local dissipation [i.e. n oc 5'(A)], and at T 00 the 
noise becomes white [u oc 5(A)]. 

For thermal equilibrium, the standard statistical mechanics result for the entropy at high 
temperature is 

5^1 + ln^ (5.1) 

k 

Obtaining this result with this formalism is a good example in its application. We will leave 
the details in Appendix B but show the numerical results in plots. Figure ^ shows a plot 
of S-Ys-z for cr = 1, /c = 1, 7o = 0.1, T = 10^. For these numbers, (|5.1|) gives S —>■ 12.513 
as 2 00, as compared with S — > 12.514 numerically at z = 100, a result indicated by the 
figure. The relaxation time, defined to be 

= 5 (5.2) 

270 

is apparent in the figure as a characteristic time over which the entropy climbs to its final 
value, while the decoherence time scale p9| 



4M7oT(t2 
is too small to be noticeable. 



^ 2.5 X 10"^ (5.3) 
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Coherent state as the state of least entropy 

We now use our entropy expression to investigate the claim that for large times the state 
of least entropy for the static oscillator is the coherent one, at least for white noise and 
local dissipation. This was shown in in the small 70 limit by using a Wigner function 
approach. 

Using our expression for the entropy S, we can plot S versus the initial squeeze parameter 
r for various times in figure ^ We have chosen k = 10, 70 = 0.1. The squeeze parameter r 
is related to a, the width of the gaussian wavefunction, by 

or. 



Note that at early times (e.g. z = 0.001), the entropy is minimised for high initial squeezing, 
as noted in p9| , fig. 1] ; this is not unreasonable since such a highly squeezed state will spread 
with time, becoming indistinguishable at later times from states which started out being 
less highly squeezed. At late times the entropy is minimised by starting with small or zero 
squeezing, i.e. an initially coherent state is the one which minimises entropy at late times. 
Thus our approach agrees with [^, and may be more useful in that it allows us to directly 
calculate the entropy at all times. 

5.1 Static inverted oscillator 

The static inverted oscillator is the simplest squeezed system. It also models the zero mode 
of the infiaton field in New Infiation ||50|- Its lagrangian is: 



L{t) = -[x^ + k^x^] (5.6) 
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Figure 2: Entropy at various times. 

Suppose this is coupled to the usual environment of harmonic oscillators in a thermal state, 
with coupling constant c(s) = 1. Then the equivalent oscillator we consider has unit mass, 
no cross term and frequency 

SO that from ( |2.11|) the sum of its Bogoliubov coefficients is (taking ti = 0) 



(5.7) 

(5.8) 
(5.9) 

r z (5.10) 

To investigate the dependence of the entropy on the various quantities in the propagator 
coefficients, we calculate these coefficients first for white noise analytically; we then calculate 
them numerically for zero temperature. 

The 6j's are independent of the temperature, and using ([A.14|) they are found to be 
(where here and elsewhere a carat will denote division by n) 



X{t) = chz — i shz 

Hence from ( |2.16D we have 

a = chz , P = —i sh. z 
so that from ( p.l4|) at late times {z oo) 



sh z 



(5.11) 



High temperature 



White noise is given by z/(s, s') = 47oT 5{s - s'), or z/(C, C') = 4:%k,^T6{C - CO; the relevant 
quantities are inserted into ( |A. 14| ) with the a^/s then becoming 



T 



an 



ai2 



2k^ sh^ z 
k'^ sh^ / 



p + e'^'/oz _ 2z - 7o ch 2z 
1 - e^^"") chz+{l + e2^«") 7o sh z 



7^g-2702 

2A;2 sh^ z 



.pg2702 - 1 + 7oe2^«' ^^2z - sh2z) 



(5.12) 



Note that 70 = 7o/«^ < 1; however if we assume small dissipation (70 ^ 1) we can write 
down large time limits of these quantities: 



an 



^70 
1-70 



ai2 



1 + 70 



^22 



^70 
1+70 



hi} 



k(±1-7o) 



-^(lT7o)2 



We can now calculate large time limits of the density matrix coefficients from ( [4 .31) : 



A 022 



B 



-bi/2 



C 



hi 



16aii 



(5.13) 



(5.14) 



These coefficients are independent of the initial conditions, which might be expected since 
the dissipation is acting to damp out any late time dependence on these initial conditions. 
So we have 

^ (5.15) 



S. 



lin 



SO that from (j^, |]T0|) 



r + 1 + In 



270T 



(5.16) 



Zero temperature 

At T = 0, the action of the environment is due to quantum effects only. If we write the noise 
in its primitive form as the usual integral over frequency then we can leave this frequency 
integration until last after the time integrations have been done. We will follow a more 
sophisticated approach in a later paper [0, but we show it here to investigate what value 
it might have. 

So we refer to ([A.14| , p.23| ), swapping the limits of integration to write 



an 



7o 



TT sh Z Jo 

7o 



duj ijj coth — - , 
IT 7o 



27r sh z Jo 



Wmax ^ ^ 

duj uj coth — - /1 1 
2T 



dC / c?C'e^°(^+^') sh {z - C) sh {z - C') coscI;(C - CO 
Jo 

(5.17) 



where 



hi = {P - c2;2 + 2e^""" + (l + 7^ + u'^^ ch 2z 

— 4e^°^ [cos Loz { chz + % shz) + 00 sin uz sh z] + 270 sh 22;| ^ 



k'' + 2u;' [1 + %) + LU 



{5.U 
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Similarly 



where 



'12 



ai2 



7oe' 



-70 z /-(i, 



TT sh Z Jo 



7- - 1 T 

da; a; cotli ^ -'12 



(5.19) 



{-2 ch^ (1 + e^^"^) - 270 sh^ (1 



+ e^"^ cos (u)^ 



— e 

2 -,2 



2702 



3 + 7^ + cZ'^ + (P - tl'^) ch 2;z] + 2cue^«^ sin cuz sh 2^} / 



and 



where 



'^22 



7oe" 



-2702 



27r sh 2; Jo 



auj uj coth — - i22 
2T 



(5.20) 
(5.21) 



'22 



{2 



=2702 



UJ 



+ (1 + 7o + cj^) ch 2z - 270 sh 2^ 

sh^]} / 



+ 4:6^°^ [cos uz {— ch 2; + 7o sh z) — u sin cl;^; 



r + 2cl^" l+7o 



UJ 



(5.22) 



With T = the coth term is set to one. Then in all cases Ojj starts at zero at ^ = 0; for low 
dissipation an, 022 quickly climb to similar constant values while ai2 climbs briefly but then 
rapidly decreases to zero. This behaviour quantitatively matches the large time limits of the 
white noise Oj^'s in ( |5.13D , even though the two calculations were done quite differently. The 
asymptotic value of an increases in even steps as we increase uOmax exponentially. So we can 
make an arbitrarily large by taking a large enough cutoff, so that it will always dominate 
D. 

In that case, with 70 ^ 1 we have at late times, using the 6j's in (|5.13| ) 



A — > a22 



B 



-bi/2 



C 



"2 



16aii 



(5.23) 



Again the coefficients are independent of the initial conditions. Since 62 is unchanged from 
the high temperature case and an, a22 tend toward constants, we now can say 



-K,e 



lin 



and so again from ( [4.7| , ^.10| ) 



S" r + 1 + In 



2^aiia22 



(5.24) 



(5.25) 



6 Scalar field in de Sitter spacetime 



We now turn to an example in cosmology, that of an inflationary universe |5^. We want to 
calculate the entropy of a massless scalar field minimally coupled to gravity in a de Sitter 
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spacetime by examining the evolution of the density matrix. As we shall see, it is a generally 
solvable squeezed system. 

Consider a scalar field $ of mass m, described by the lagrangian density 



C 



'9 



(6.1) 



coupled by ^ to the curvature R = Q{a? /a^ + a/ a) of a spatially flat FRW universe with 
metric 

ds^ = dt" - a\t)Y,{dx'f (6.2) 

i 

In conformal time r] = J dt/a, the conformally- related field x = a$ is described by a La- 
grangian density 



C 



X' 



E2 o ' 2(22 " I 
X,i -2-XX -X [ma ^ + 6^ 



/2 



.a 
a 



(6.3) 



Decomposing the field into normal modes k with amplitudes qk, the Lagrangian can be 
expressed as 



,/2 



'2 



2-qq' - gM + m^a' - + 6^— 



2„2 



(6.4) 



Inside the brackets if we add a surface term of 6^(g^a'/a)' to eliminate the a" term (for 
justification, see [0]), we get a new lagrangian: 



Lnewiv) — E 2 



+ 2(6^ - l)-qq' - q' ( k' + m'a' + (6^ - 1)^' 



a 



a" 



(6.5) 



For a massless minimally coupled scalar field in de Sitter space. 



Lnewiv) — E 



/2 I ^ / 2/72 

q + -qq - q [k 

V V 



(6.6) 



We also use a spectral density of the form 

I{u;,r],r]') 



27o u 



(6.7) 



7tH y^rjrf 

so that c{ri) = l/y/^Hrj. This form of spectral density will be justified in a later paper [H5 



although for now we note that it does not make the equation of motion for X any harder 
to solve than if we had used a static coupling. Since '~fo/H is dimensionless we rewrite it as 
c [not to be confused with c{ri)]. Incorporating the bath gives the equivalent oscillator with 
M = 1,£ = l/r] and frequency, from ( |A.6[ ), 
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Also we choose k = k to simplify the equation of motion. With z = kr] we can write this 
together with its initial conditions from (|2.9| , |2.11| , |2.12|) as 



X"{z) + 
X(z,) = 1 



2 + c2 



X = 



X'izi 



- l/zi (6.9) 

where z < 0. The solution of this equation can be constructed using Bessel functions whose 
index is a function of c; however since we are interested in small c we take the solution to 
be approximately that of the same equation but with c set to zero. This simplifies things 
greatly: 



X(z) 



where 



2zi 



1 



J{zi-z) 



(6.10) 



(6.11) 



We can further simplify X by using a very early initial time, setting Zi —>■ — oo. We also 
disregard the phase in the resulting expression for X, since this is not expected to make any 
difference to physical quantities. In this case we obtain a new function which we rename X: 



X{z) 



1 - - 1 e-'' 
z, 



The Bogoliubov coefficients can now be found from ( 2.16| ): 



and so from ( |2.14| ) at late times 



In \z\ 



(6.12) 



(6.13) 



(6.14) 



This result was also obtained in using a different formalism. 

First we calculate the fej's. Since we are only interested in late times we can work to 
leading order in z (although with hindsight we include some next higher order terms which 
will be needed later). Using (|A.14 ) we find 

ck/z + kz + 0{z^) 
Tk\z\ 



hi 



hi) 



\l^c\ \±c 

^7, 



64 = {c+l)k/zi + kz'^ /?, + 0{z^) 
and for the aj/s we need the following expressions, calculated from ( |6.12| ): 

Im [X{z)X*{C)] (1 - z/0 cos(C -z)-{z + l/Q sin(C - z) 



(6.15) 



exp 7o 



Im X{z) 

Im [X(C)] 
Im X{z) 



cos z + z sin z 



^ + sinC 



M 



-dC 



+ sin z 
; exp 



-7o 



M 



(6.16) 
(6.17) 
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6.1 Entropy 

High temperature 

We begin by writing 



V = 4cc\s)T5{s - s') 
-Ack'^T 



(6.18) 



We calculate an here and leave the details of ai2, 022 to appendix 0. First, ( |A.14| ) gives 



ail 



2k 



2cT / dC 



^Zi J Im X{z) 

Vim [Xiz)X*iC)]) 



Zi 



Im Xiz] 



1 



(6.19) 



We wish to investigate the dependence of the a^j's on ^ as z ^ 0, and so we now separate 
each integral into a sum of two parts. The first is gotten by integrating in to some constant 
A close to 2;, while the second integral contains the z upper limit: 



an = 2cT 



+ 



'^V lm [X(^)X*(C)] y 1 
I Im X(z^ ] 



< 



(6.20) 



It's only necessary to work to leading order in z. We need the following expressions: when 
only z ^ we have the z dependence in the integrands as 



Im \X{z)X\0\ 
Im X{z) 
Im [X(C)] 



Im X{z) 

while if both then to leading order 

Im \X{z)X\Q\ 



cosC - sinC/C + 0(^2) ^ /i(C) + 0(^2) 
^(cosC/C + sinC)=;^/2(C) 



Im X{z) 

We are now in a position to write 



(-C + ^70/3 



Im [X(C)] 

Im xr^) 



vc 



(6.21) 



(6.22) 



On oc cT 



cT(0(1) + 0| 
cTO(l) 



-2c+5 



(6.23) 



since we have taken c to be small. A similar approach gives the following results for 012, 022 
(details can be found in appendix |C|): 



ai2 = cT 0\z 



c+l 



022 = cT 0(1) 



(6.24) 
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Since T is large, an dominates D while 022 dominates A\ so we have 



A 022 



B 



-h/2 , C 



16aii 



(6.25) 



These of course have the same form as for the static oscillator case, although it's by no means 
clear whether such a fact could have been deduced from the general expressions for the a^j's. 
We now have 

Slra^ , ~^ = 0|Z|1-^ (6.26) 



and using (^]^, |6.14| ) we can write 



5 — > (1 — c)r + constant 



(6.27) 



Finite temperature 



Here we leave the frequency integration until last as was done for the static oscillator. The 
integrals can then be done in the same way as in the last section, although some subtleties 
are present in this case (see appendix 0). We finally obtain 



an = ck 0(1) , ai2 = ck 0\z^l'^ , 022 = ck 0{z) 
Again since we integrate over cj, an will be large and so dominate leading to 

A-a22-^ , 5 ---61/2 , C 

and so 



"2 

16aii 



Then with ( [4.7| , |6.14| ) we have 



S ^ (1/2 — c)r + constant 



(6.28) 

(6.29) 
(6.30) 
(6.31) 



7 Discussion 

In the last two sections we calculated the entropy of two physical and exactly solvable 
squeezed systems; an inverted harmonic oscillator and a scalar field mode evolving in a de 
Sitter inflationary universe. Our aim was to compare these results, based on our rigorous 
quantum open system framework, with that of the previous more ad hoc approaches de- 
scribed in the introduction. We must bear in mind that these previous results referred to a 
field mode that could be split into 2 independent sine and cosine (standing wave) compo- 
nents. We should therefore expect a result of 5* = r (rather than 2r) if we are to agree with 
previous work. 

For the inverted oscillator, in both temperature regimes with low coupling, we obtained 
S — >■ r+constant. In the de Sitter case, the high temperature result is S* ^ (1— c)r+constant. 
These three examples certainly do confirm the ad hoc approaches to calculating entropy 
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that have been used by others. However at lower temperatures the de Sitter entropy is 
S (1/2 — c)r + constant. This last result requires us to look more closely at A and C 
which together give the entropy. 

From ( [4.6| , [4.7|) , and neglecting the added constants which are always implied, we find 



that in the high squeezing limit the entropy behaves as 

S ^ -hiA- -InC. 
2 2 

When the system-environment coupling is small, all of the above cases give —1/2 In C r, 
which is the expected result. The dominant contribution to C always comes from 62 in the 
high squeezing limit. This parameter is determined by the squeezing of the system and is 
essentially independent on the nature of the environment and its coupling to the system. We 
can therefore conclude that the In C contribution to the entropy represents entropy intrinsic 
to the squeezed system itself. This is in agreement with the previous results and should 
also be true quite generally for squeezed systems. However these results cannot but fail to 
take into account the contributions to the entropy from the In A term. This contribution is 
determined by the aij factors which strongly depend on the nature of the environment and 
its coupling to the system. There is, a priori, no reason to expect this contribution to be 
small, a point illustrated by our finite temperature de Sitter example for which we found 
1/2 InA — >■ — r/2. This highlights the danger in using the previous ad hoc approaches to 
entropy of squeezed systems. The critical point is that the entropy of a system depends not 
only on the system itself but also on the nature of the environment it is coupled to. 

In conclusion, approaching the problem of entropy and uncertainty from the open system 
viewpoint as we have demonstrated improves on the earlier work in that it makes explicit 
how their dependence on the coarse-graining of the environment and the system-environment 
couplings. It also clarifies the relation between quantum and classical descriptions - it is 



through decoherence that the quantum field becomes classical p2| , |52| . These issues are im- 
portant as they rest at the foundation of statistical and quantum mechanics. (For a discussion 
of the deeper meaning of the dependence of persistent structures on coarse-graining, see ||18[|.) 
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A Details for section 



A.l Calculating ui 



V2 



Now we are in a position to solve ( p.l7| , |3.18| ) for ui V2- First consider ( p.l7| ). We treat 
the integral of a delta function and its derivative in the following way: use a smooth step 
function (i.e. 9{0) = 1/2) to write (xi > xq) 



Xl 



Xo 



f{x)6{x — a) dx 
f{x)6'{x — a) dx 



f{a) 6{xi — a) 6{a — a;o) 
-/'(a) 9{xi — a) 9{a — xq) 



(A.l) 
(A.2) 



These relations can easily be proved by checking the five cases individually, of a < Xq, a = xq, 
Xq < a < Xl etc. Note that treating the delta function in this 'smoothed' way eliminates 
the need for the frequency renormalisation in |^9[. This smoothing essentially just defines 
6{x)dx = 1/2 (see e.g. ^Tj for a discussion of this). 

Hence ( p.l7| ) together with ( p.23| ) becomes (with u being either ui or U2) 



'M 27oc2 



Now define u by 

in which case it follows that 



u = u exp 



70 



MS 



. Mis'] 



+ £ 



ds' 



M , 



u = 



M \ M M2 



u = 



(A.3) 



(A.4) 



(A.5) 



Comparing with ( p.ll| ), we recognise this as just the equation of motion of an oscillator 
with mass M, cross term S and an effective frequency 



2 



M2 



(A.6) 



So, we are in a position to describe our system in terms of an equivalent system. Hence 
we know a solution for uls) — it is the sum X of the Bogoliubov coefficients for this new 
system. So we write (with (71,(72 constants to be determined) 



u{s) = exp 



-70 



—ds' 

. M 



[(7iX(s) + (72X*(s)] 



(A.7) 



By including the boundary conditions for ui and U2 we obtain 



ui(s) 



U2(S 



exp 



exp 



-70 



7o 



M 



M 
ds' 



:ds' 



Im [X{t)X*{s)] 
Im X{t) 

Im X{s) 
Im X{t) 



(A.8) 



26 



This tying in of the propagator formahsm to the language of squeezed states (such as Bogoh- 
ubov coefficients) will be very useful for relating the entropy of a field mode to its squeeze 
parameter r. 

In the same way that we solved ( |3.17|) , eqn ( p.l^ ) becomes 



(M 27ocn . ^(r,2^M£ ■ 27oCc\ 
v{s) + 77 r-;— \v+\il + ^rr + £ tt- I f = 



M 

Now write 

V = V exp 
and just as for the case of u we have 



M 



M 



V H V + H V£ - U = 

M V M M2 / 



(A.9) 



(A.IO) 



(A.ii: 



So now Vi and V2 can also be written as combinations of X and X* . Including the boundary 
conditions we eventually obtain 



= exp 
V2{s) = exp 



-70 / j^ds' 



Im [X{t)X*{s) 
Im X{t) 

Im X(s) 



Im X{t) 



(A.12) 



A. 2 Calculating an ^ 64 

To facilitate our calculations we introduce dimensionless parameters for time 

Z = Ht , ( = KS 

X{z) = X{t) etc. 



(A.13) 



and a carat will denote division by k, e.g. 70 = 70//^- Note that t is the lagrangian time, 
which isn't necessarily cosmic. 

Now we are able to calculate the propagator. Making use of ( p.l6| , |3.15| ) we obtain 



2k2 



dC / dC 



, ^,0/,^ #.C" Im [X{z)X*{C)] ^(^^ ^70 4' g.c" Im[X(^)X*(CO] 



Im X(z 



Im Xf^) 



K J Zi J Zi im (^zj im (^zj 

1 /"^ ,> r --Yn f'^dC'Im X(C) £.dc"lm. X(C') 



bi{z,Zi) 



-70 «:c2(^) + «:M(z)— — + M(;.)£:(z) 
Im A u 
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Im Xiz) 



-70 KC [Zi) + K 



Re X{z) 
Im Xf;z) 



+ M(zi)£:(z, 



(A. 14) 



B Entropy of a static oscillator in a thermal bath 

The lagrangian for the static oscillator with unit mass is given by 

L = 1 _ ^2^2^ 

From ( |A.6|) with M = c = 1 the effective frequency is 

^eif = — 7o = 

Then the equation of motion for X is, from (|2.11|) with Vt ^ VL(.jj 

X + K^X = 
X(0) = 1 , X(0) = -ZK 



which leads to 
with z = Kt. Then 



X{z) = e-'' 

lm[X{z)X*{C)] sin(;z-C) Im X(C)] sinC 



Im X{z) 



sin z ' Im X(z) sin z 



e •'^i ' = e 
with noise for T ^ 00 being white: 

KC, = 4«:7oT5(C - C) 



Then an ^ 64 follow: 
an 

ai2 

022 



T e^^o^ - l-7o sin 2^-72(1- cos 2x 



sin^ z 



2(1 + 70^) 

2T — C0S2; sh7oZ + 7o sin z ch7oZ 



sin^ z 



1 + 70^ 

T -e-2^"^ + l-7osin2z + 72(l-cos22) 



sin^ 2; 



2(1 +7o') 



= /t(-7o ± cot 2;) 



sm z 



(B.l) 

(B.2) 

(B.3) 
(B.4) 

(B.5) 

(B.6) 
(B.7) 

(B.8) 



(B.9) 
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To evaluate S, we need A and C; in turn for these we need an 
from (10^) .1 

The oscillator is assumed to be initially in its ground state 

^2 



ip{x, 0) oc exp 



-X 



4^2 



so that its density matrix is 



and in (14.11) we have 



p{xx' 0) oc exp 
1 



— — 



4^2 



4a2 



The reduced density matrix evolves into 



, X = 
with 



64. These are calculated 

(B.IO) 

(B.ll) 
(B.12) 



A 
B 

C 
D 



0.22 + 



D 



1 



8(t2 



an 



&3 + «12&3&4 



"12 
2a2 



2 2L) ^ 

ft2 
"2 



ai2l 



8Da^ 



1 , 2an , 
H ^ + 0. 



It's by no means trivial to show that the entropy calculated using these expressions does 
indeed tend toward ( ^.1|) , and in particular the esc 2; terms in the a^j's and fej's mean their 
values can diverge depending on the time. But this divergence cancels out when physical 
quantities are measured, as we can see by verifying numerically that our entropy really does 
tend toward the usual asymptotic value at late times (Fig. 1). 



C Calculation of a^j's in section |6| 
de Sitter with high temperature 

Here we evaluate the a^j's leading to ( |6.24| ). We are using the following small z,( approxi- 
mations: 



Im [Xiz)X*iC)] 
Im X{z) 

Im X(C) 
Im X{z) 



^ cosC -sinC/C + 0(^2) = /i(C) + 0(^2 



z(cosC/C + Sin = ^/2(0 

z/C (C.I) 



■^Various notations exist describing these results; see for example jlj, ^ To compare with ||2^, 

eqn 2.2.7] is a matter of carefully transcribing the notation; key things to note are that X = I], Y = —A; 
here we have taken xq = Pq = 0; ||2^, eqn 2.2.6c] should have an an in place of the 022; the bi's in |2^ are 
written explicitly in if^] via [pi 3.11]; p|, 012] equals p^, ai2 + 021]; p|, 70] equals p^, 7o/2]. 
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Firstly, 

K J Zi J Zi \ Zi 



Im X{z) 



-c 



.C7 ImX(z) 



4cr 



C Im [X(z)X*(C)] 1 /zy imX(C) 



oc cT\z\ 



f'dc \C\-''-'MOzMC)+ fdQ \C\-"'-\-C + z'IQ 

J Zi J X 



z 

3C 



c+l 



provided c < 1/2. Finally, 

= w^Jz/^Jz/^ [c) i^xR^^^^-^Hc 

= cro(i) 



Im X(^) 



Uv' llmXC^); -C 

v A 



(C.2) 



(C.3) 



de Sitter with finite temperature 

We leave the frequency integration until last: 

2c 1 



2c A;^ i-oo 



= / a; coth — - cosa;(s — s') 
> Jo 2T ^ ' 



x/CC 



= duj Cj coth — - cosa)(C — C') 
' Jo 2T vs s y 



The ajj's are 



On 



* TT io 



z;'' — / du uj coth — X 
2T 



(C.4) 



^^c r< ici-v- ^- cos.(c-co icr-v^ ^- f (-)f (CO] 

;i Jzi Im A (zj Im A (zj 



=/ll 



ai2 



[ZiZ) / du! cu coth — X 

TT Jo 2T 



r''ci:<icr--5=^i™ico..(c-c')ic'i 



-c-l/2 



Im X(C') 
Im X(z) 



022 



ck 



TT JO 



2; — / da; a) coth— X 
2T 



=722 
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Using the expressions from ( |C.1| , |6.17| ) the first of the inner integrals becomes 



'11 



J Zi 



+ 
+ 



dc coscocicr^~'/'(-c"+^vc')/3 



-c-l/2/_^2 



(-C + ^7C)/3 



f^dC COSiuC IC 
Jzi 

+ fdc coscoic-c) \cr-'^\-c"+zVcm 

J A 



n-c-l/2 



/i(C') 



(C.6) 



We now have a difficulty. In order to get a reasonably useful analytic result, it will be an 
advantage to replace the cos uj{C,—C,') term in the fourth integral above by something simpler. 
We will have competition between uj increasing in the frequency integral versus z decreasing 
in time. Suppose then we use a frequency cutoff ujmax- In that case we can approximate 
cosa)(C — C') foi' C; C ~ by choosing Umax such that coscj(^ — ^ 1 in the fourth integral. 
This will be true provided 

UJma. < -1/A (C.7) 

However now we don't expect our result to necessarily agree with the high T result found 
in ( |6.23| ), since there we had taken uJmax — ^ cc, which was made possible by the use of the 
delta function. 



At this point we refer to the discussion of the high temperature limit in There it is 
shown that the high temperature (delta function) regime is that for which uJmax T and 
^max ~^ cxD. This absence of a cutoff in the high temperature limit is usually not stressed, 
but it forms the most relevant fact here. In general we must impose a cutoff for all finite T 
values, otherwise the frequency integral is not well defined — unless T — >• oo. So we conclude 
that the regime for which our analysis is valid here is T < Umax- 

With the last cosine set equal to 1 as before, these integrals are all 0(1) and therefore so 
is ail. Next: 



'12 



dC COSLuiC-C) IC 



/|-c-l/2 



f2iC)z 



+ 



dC cosc2;C IC 



+ / dc\cr-'/\-e + zyo/3 



[ dc cosuc\cr-'^'f2{c')z 

J Zi 



+ / < coscI;(C-CO IC 



/|-c-l/2 



Evaluating these integrals gives /12 = 0\z\ ^+1/2 go that ai2 = 0\z\^/'^. Lastly, 



(C. 



'22 



^cicr^-'/V2(c)^ 



dC cosuj{C-0 IC 



-c-1/2 



/2(C')^ 



+ 



dC cosuC IC 



/|-c-l/2 
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J \ J Zi 

+ rdc COS coic-c) \cr~'^'z/c' 

= 0|z|-2<^+^ (C.9) 
so that = 0{z). 
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